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Abstract 

The ratio between the anomalous magnetic moments of proton and 
neutron has recently been suggested to be connected to the ratio of proton 
momentum fractions carried by valence quarks. This ratio is evaluated us- 
ing different constituent quark models, starting from the CQM density dis- 
tributions and calculating the next-to leading order distributions. We show 
that this momentum fraction ratio is a sensitive test for SU(6)-breaking 
effects and is a useful observable to distinguish among different CQMs. We 
investigate also the possibility of getting constraints on the formulation of 
quark structure models. 
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1 Introduction 



The static properties of baryons are an important testing ground for QCD based 
calculations in the confinement region. However, different CQMs[|l], ||, ||, ||, ||, 
are able to obtain a comparable good description of the low energy data, so 
that it is difficult to discriminate among them. A fundamental aspect of the 
theoretical description is the introduction of terms in the quark Hamiltonian 
which violate the underlying S't/(6)— symmetry. It is therefore important to find 
out observables which are sensitive to the various SU(6)-breaking mechanisms. 

In this respect, the relation proposed recently by Goeke, Polyakov and Van- 
derhaeghen between the anomalous magnetic moments of the proton and the 
neutron and the proton momentum fractions carried by valence quarks, Mf""', 
might be a good candidate for testing SU(6)-breaking effects and can lead to 
important constraints on the models for the structure of the nucleon. 

Quark models are able to reproduce in a extraordinary way the static low en- 
ergy properties of baryons with very few parameters and this gives us confidence 
that they are a good effective representation of the low energy strong interaction 
dynamics. The QCD based parton model reproduces in a beautiful way the 
dependence of the high energy properties even with naive input. However the per- 
turbative approach to QCD does not provide absolute values of the observables; 
one can only relate data at different momentum scales. The description based 
on the Operator Product Expansion (OPE) and the QCD evolution require the 
input of non-perturbative matrix elements which have to be predetermined 
and therefore the parton distributions are usually obtained in a phenomenological 
way from fits to deep inelastic lepton nucleon scattering and Drell-Yan processes. 
The basic steps are to find a parametrization |10| which is appropriate at a suf- 
ficiently large momentum Qq^? where it is expected that perturbation theory is 
applicable, and then QCD evolution techniques are used in order to obtain the 
parton distribution at higher Q^. Using these parametrizations a large body of 
data is reasonably described, even if at the origin this parametrization is purely 
phenomenological. 

Cluck, Reya and Vogt |Tl| started from a parametrized distribution of partons 
at a very low scale /z^q, which resembles that of a naive Quark Model of hadron 
structure, in the sense that the contribution of the valence quarks to the structure 
function is dominant. As suggested by Parisi and Petronzio |]r2| , the hadronic 
/Xq scale is defined such that the fraction of the total momentum carried by the 
valence quarks is unity. This procedure opens the possibility of using Constituent 
Quark Models as input in order to calculate the nonperturbative (twist-two) 
nucleon matrix elements, as proposed by Jaffe and Ross E 



The scheme developed by Traini et al.|14| takes into account all these aspects: 
it uses as input the quark model results in order to determine the non perturba- 
tive matrix elements at the hadronic scale [|12], then an upwards NLO evolution 
procedure at high momentum transfer (Q^ = 10 GeV^) is performed . 



2 



Starting from three different Constituent Quark Models [|T], we have 

calculated the parton distributions at the hadronic scale and we have evaluated 
the ratio of the proton momentum fractions carried by valence quarks. A NLO 
evolution has been performed up to = 10 GeV^. 

All models give a good description of the spectrum and have been used also 
to describe various observables (elastic and inelastic form factors, strong decays). 
In particular, the different results for the electromagnetic transition form factors 
indicate that the models have a quite different Q^-behaviour. However, as we 
shall see later, the ratio of the proton momentum fractions carried by valence 
quarks is independent of the scale Q^, therefore we expect that the study of this 
relation will give important information on general aspect of CQM. 

The paper is organized as follow. In Section 2 we review in a critical way the 
new relation as found in Ref. § between the ratio of the anomalous magnetic 
moments of the proton and the neutron and the ratio of the proton momentum 
fractions Ml""'. In Section 3 the unpolarized parton distributions are evaluated, 
at the hadronic scale, using different CQMs, and an evolution procedure is per- 
formed and then in Sect. 4 the ratio of the proton momentum fractions carried by 
valence quarks is calculated as a function of and compared with experimental 
values and with the results of the models for the ratio of the anomalous magnetic 
moments. 



2 Ratio of proton momentum fractions carried 
by valence quarks 

In Ref. a relation has been proposed between the ratio of the proton and 
neutron anomalous magnetic moments and the momentum fractions carried by 
valence u- and d-quark distributions, as follows : 

kP _ 14 M2 + Ms"""' 
ft" ~ 2 Ma''""' + M^"""' ' 

with the proton momentum fraction carried by the valence quarks defined as 

Mf"' = 1^ dxxq^ai{x). (2) 



In Fig. ^ we show the scale dependence of the rhs of Eq. (|lD, which we shall hence- 
forth denote with R, for various recent parametrizations of next-to-leading order 
(NLO) and next-to-next-to-leading order (NNLO) parton distributions. Fig. |l] 
shows that the scale dependence drops out of the rhs of Eq. (|I]), although the 
numerator and denominator separately clearly have a scale dependence. Fur- 
thermore, it is seen from Fig. |I], for all NLO and one NNLO parametrizations 
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of parton distributions, that the relation of Eq. (|1]) is numerically verified to an 
accuracy at the one percent level! In particular, the most recent MRSTOl NLO 
frg, the MRSTOl NNLO 0, and the CTEQ6M NLO |8l parton distributions 
(which appeared after the writing of Ref. nicely confirm the finding of Ref. ^j. 
Although the relation Eq. (|1|) was originally derived within a parametrization of 
generalized parton distributions, it is in fact completely independent of such a 
parametrization, as the rhs of Eq. (|l]) is expressed in terms of moments of forward 
valence quark distributions alone. 

The above observations from phenomenology suggest that Eq. (|l]) holds and 
that the unpolarized valence u— and (i-quark forward distributions contain a 
non-trivial information about the anomalous magnetic moments of the proton 
and neutron. It is the aim of the present work to investigate the relation of 
Eq. (H) in different quark models. 

Let us firstly consider the simplest quark model, with exact SU (6) symmetry. 
In this limit, M^^""' = 2 Ai^^""', and = - = 2, so that one immediately verifies 
that Eq. (|I]) holds. 

In reality, the ratio of anomalous magnetic moments deviates from the SU{6) 
limit by about 6.5 %. The smallness of this deviation is the main reason why 
constituent quark models are quite successful in predicting nucleon (and more 
generally baryon octet) magnetic moments. In quark model language, the rela- 
tion of Eq. ([1|) implies that the small breaking of the SU (6) symmetry follows 
some rule which is encoded in the valence quark distributions. In particular, it is 
interesting to investigate a possible correlation between the ratio of valence d— 
and M-quark distributions, and the ratio of proton to neutron anomalous mag- 
netic moments in different models. To this end, we turn in the next section to the 
calculation of parton distributions in quark models with different SU{6) breaking 
mechanisms. 



3 Parton distributions from quark models 

The approach, recently developed by M. Traini et al. for the unpolarized distri- 
butions ||l4|, connects the model wave functions and the parton distributions at 
the input hadronic scale through the quark momentum density distribution. In 
the unpolarized case one can write the parton distributions ||14|| : 

= / d'k n,(|k|) - ^) (3) 

where fc+ is the light-cone momentum of the struck parton, and nq(|k|) represents 
the density momentum distribution of the valence quark of q-flavour: 

n,(|k|) = {N, J, = +1/2I ^ 6{k - k,)\N, J, = +1/2) 

1=1 ^ 
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nrf(|k|) = (iV,J, = +l/2|^ 



5(k-k,)|iV,J, = +l/2) 



(4) 



is the third component of the isospin Pauh matrices, ki is the momentum of 
the ith constituent quark in the CM frame of the nucleon, \N, = +1/2) is the 
nucleon wave function (in momentum space) with = +1/2 component. 
Using k+ = ko + kz, one can integrate eq. 13 over the angular variables and get: 



gy(x,/io) 



27rM 



d|k||k| r2,(|k|), 



(5) 



where 



M 
2 



X 



1 — X 



X 



1 - X VM, 

M and are the nucleon and (constituent) quark masses respectively. 
Eq. (H) can be applied to a large class of quark models and satisfies some impor- 
tant requirements: it vanishes outside the support region < x < 1 and it has 
the correct integral property in order to preserve the number normalization. 

In the present section we shortly illustrate the evolution procedure we have 
been using. Even if alternative factorization schemes have been investigated, we 



remain within the MS renormalization and DIS factorization scheme(see [|T5l and 
references therein). In this case the moments of the F2 proton (neutron) structure 
functions have the simple expression 



Yl ^?(^ ^(^' Q"^) + ^ ^(^' <5^))n 



q=u,d,s 

2 



+ (-)^ {xq3iQ'))n + ^ {xqsiQ'))n + I {xnQ'))n 



(6) 



IS a 



where + (— ) refers to proton and neutron respectively; S = Y,q{(l + o) 
singlet component and = u + u — {d + d), gg = u + u + d + d — 2{s + s) are 
nonsinglet (NS) contributions. The Wilson coefficients C^^''^ and C^^'^ , in the 
MS renormalization and factorization scheme, can be found,e.g., in Refs. P^, 

The NLO evolution of the unpolarized distributions is performed following 
the solution of the renormalization group equation in terms of moments, i.e. 
(/(Q^))n = Jo dx f{x,Q'^) x"^^^. Since, in our case, the starting point for the 
evolution (/Xq) is rather low, the form of the equations must guarantee complete 
symmetry for the evolution from Hq to ^ /Xq and back avoiding additional 
approximations associated with Taylor expansions and not with the genuine per- 
turbative QCD expansion |jT5|]. In particular for the Non-Singlet sector we have 



{qNs{Q'^))n 










2/3o 






l,n 

7jvs 






2/3o 




47r 



(7) 
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where Jn'P'^ ci^^ the anomalous dimensions at LO and NLO in the DIS scheme 
Q, and Po, Pi the expansion coefficients (up to NLO) of the function 
/?o = 11 - 2/3 Nf, pi = 102 - 38/3 Nf for Nf active flavors. Eq. (0) reduces to 
the more familiar form (e.g. Ref . ||TT|, |23| ) 



{qNSm)n 




l,n 

Ins 



2Pi , 



4n 



after performing a Taylor expansion for both 



47r 



•C 1 and 



< 1. 



The A's values are suggested by the analysis of Gliick et a/.[TT|, Q;s(/io)| 



I NLO 



is obtained evolving back the valence distribution as previously mentioned, and 
/Xq is found by solving numerically the NLO transcendental equation 



In 



^*-NLO 



Poot 



47r Pi 

Pi 



Poa 







(9) 



which assumes the more familiar expression 

«s(Q') 1 



An 



/ Po Inln(QVA^Lo) 

/3oln(QVALo) V Pi In(QVALo) 



(10) 



only in the limit ^ A^lq; (an interesting discussion on the effects of the 



approximation (|T^) can be found in Ref.j 

The actual value of /Xq is fixed evolving back (at the appropriate pertur- 
bative order) unpolarized data ffis, until the valence distribution x V{x, Hq) = 
X'Uy(x,/ig) + X dv{x, nl) matches the required momentum (/ dx xV{x, fi^) = 1). 
The resulting NLO (LO) parameters are 



a 



in 



NLO 



LO 



0.142^ 
0.290. 



NLO 



LO 



0.094 GeV% Anlo 
0.079 GeV^ Alo 



248 MeV ; 

232 MeV. (11) 



We discuss the results obtained using different models for the valence quark 
contributions, namely the Isgur-Karl (IK) model |jl| , which has been largely used 
in the past to study the low-energy properties of hadrons and also deep inelastic 



polarized and unpolarized scattering|15], a hypercentral Coulomb-like plus linear 



confinement potential model inspired by lattice QCD and an algebraic 
"'^The 7]y^ are redefined in the DIS sclieme in such a way that the Eq. holds, i.e. 7^^^ — > 
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model the wave functions of the last two models give a rather good descrip- 
tion of the electromagnetic elastic and transition form factors IQ p5| pT]] . 



1) The well known Isgur Karl model is based on a harmonic oscillator potential 
plus a One-Gluon-Exchange-hyperfine interaction which is responsible for the 
SU (6) breaking of the symmetry. The nucleon wave function is written as a 
superposition of SU (6) configurations, that is 

\N) = as\5Q,0+) + a's\5Q',0+) + aAf|70,0+) + a^|70,2+). (12) 

In particular we discuss the result for the Isgur Karl model(IK) as = 0.931, a'g = 
—0.274, ttM = —0.233, = —0.067 and also for a simplified model where only the 
as and qm (or qd) coefficients do not vanish. The contributions from the SU(6) 
breaking components come from the amplitudes a's, au and an of the |56',0+) 
|70,0"'") and |70,2+) multiplets, since without the OGE-hyperfine interaction 
ci's = ci][,f = ci£) = 0. 

The corresponding momentum density distributions are 



^n(|k|) 




3/2 



+ a 



/ 2 



3 k^ 3 k^ 

2^ + 4^ 



+ a 



k^ ^k^^ 



- asa'gA\/3 — - 1 + asaMVQ 




2a^ 




(13) 



nd{\k\) 



with 



and 



+ 



.gs 3/2 

2 a%3/2 y 2 j 



^5 



I / 



s 



3k2 3k^^ 
2 a2 + 4 ^4 



I 2 



Ik^ 3k^^ 

4 a2 ^ 8 




asa's2Vs ( ^ - 1 ) - asttM^Q ( ^ - 1 



+ 



asttM 



.1V2 




J nu{\k\)dk = 2 J ndi\k\)dk = 2 



n(|k|) = n„(|k|)+n,(|k|) . 

The ensuing M2 momenta for the u and d quarks in the proton are reported 
in Fig. 2. The scale dependence for the single momenta is quite smooth apart 
from the low Q^-values. 



+ 



(14) 
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2) The hypercentral Constituent Quark Model (hCQM) is based on a Coulomb- 
like potential plus a linear confining potential to which a OGE-hyperfine inter- 
action is added. The difference with the IK model is mainly in the spatial wave 
functions which are not gaussians, but are more spread out and are obtained by 
numerical solution of the 3— quark wave equation. Moreover, the nucleon state is 
written as a superposition of five SU(6)- configurations: 

\N) = a5|56,0+) + a5|56',0+) + a'^|56",0+) + aM|70,0+) + 0^170,2+) (15) 

with as = 0.9997, a's = 0.0217, = 0.0041, om = 0.0038, = -0.0012. With- 
out the 0GE-SU(6) breaking term a'g = a'g = au = = 0. The resulting 
momentum density distribution contains higher momentum components in com- 
parison with the h.o one. 

2_bis) The SU(6) invariant hamiltonian is left unchanged, while the SU(6)- 
breaking mechanism is provided by a spin- and isospin-dependent interaction ||26|| . 
Here also the content of high momentum component is greater than in the h.o. 
case. 



3) In the model proposed by lachello et al.p| the hamiltonian consists of a 
part corresponding to the vibration and rotation of a top to which a Giirsey- 
Radicati spin and isospin dependent term is added. The Giirsey-Radicati term 
is diagonal with respect to the SU(6)-configurations, so it splits but does not 
mix the ^[/(G)— configurations. An S't/(6)-breaking mechanism is implemented 
in a phenomenological way considering different u and d charge distributions 
27| , which correspond to different u and d effective charge radii. In this way the 



nucleon elastic form factors are obtained folding the top form factors with the u 
and d charge distributions (assumed to be exponential-like) and the results have 
a dipole behaviour. 

Also in this case the momentum density distribution contains high momentum 
components and one can imagine that this will strongly influence the results. 
The unpolarized valence quark distributions are given by 

8 Kf 
^ ^1 + k 



where: 

Mu = 2 J\fd = l 



a?} = 0.258 fm = 0.285 fm. 
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The validity of Eq. for the model 2_bis is analyzed in Fig. 3. The two 
members are equal within 0.2 %, although the K-ratio differs by about 7 % from 
the experimental value (~ —0.937). 

Similar results, reported in Table I, hold for the other models, with the excep- 
tion of the U(7) model, where the k- value is correctly reproduced by construction, 
while the equation is violated up to a few percent. 

In order to test if this feature depends on the choice of the CQMs or is a 
general characteristic, we have used the analytic expression supplied by the Isgur- 
Karl model and tried to reproduce the experimental value of the two ratios by 
leaving the amplitudes a'g,aM and ao free. One can also vary the h.o. constant 
a, with being a measure of the confinement radius. The Q^-behaviour of 
the I.K. model is unrealistic because of the gauss-factors, however also in this 
case the ratio is quite scale independent. The procedure of fitting the amplitudes 
corresponds to introduce implicitly quite different hamiltonians. The anomalous 
magnetic moments have the following expressions: 



2(1- 
-2fl 



I'M) - 



4a 



D 



3/2 



(17) 



If one adopts a model where the only SU(6) breaking comes from the om, it 
is immediately seen from equation (|T^ that the K-ratio is exactely equal to -1, 
like in the SU(6) limit. The crucial quantity seems then to be the a/) amplitude. 
Assuming that the D-wave amplitude is the only SU(6)-breaking term {D-model), 
we have that: 



2q| - al 
-2a| - 1/2 al 



-0.937 



if = 0.955 and = 0.295. Calculating the rhs of Eq. (|I]), which we refer as 
R in the following, with these two values of the parameter and varying a in a 
quite large interval, the best value obtainable is i? = 0.9988, with a = 2.1 fm~^, 
differing by about 7% from the K-ratio. Finally, leaving completely free the 
amplitudes a'g, qm and a^) in order to fit the /t-ratio and R separately, the resulting 
amplitudes turn out to be complex. 

Therefore, the proposed Equation (|l|) seems to be valid (up to few percent) 
for all Constituent Quark Models provided that the SU (6)- violation is not too 
strong, but both values are quite far from the experimental value of the K-ratio 
of —0.937. If one tries to force the SU(6)-violation to reproduce the experimental 
value, one is apparently faced with too strong constraints coming from the CQM 
itself. This is a possible indication that the degrees of freedom introduced in the 
current CQM may be inadequate since one has to take into account pion cloud 
effects 
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4 Discussion and Conclusions 



The relation Eq. (|1|) between the ratio of the proton and neutron anomalous 
magnetic moments and the momentum fractions carried by valence quarks, Mf""' , 
is exactly verified in the SU(6)-invariant limit, where both are equal to -1. 

In the currently used Constituent Quark Models, SU(6) violations are intro- 
duced in different ways (One-Gluon- Exchange interaction, spin and/or isospin 
dependent terms, Giirsey-Radicati mass formula, One-Boson- Exchange ...). Such 
SU(6) violation is necessary in order to bring the anomalous proton and neutron 
magnetic moments closer to the experimental values or to reproduce important 
features of the spectrum, such as the N-A mass difference. 

In all the models we have considered in this paper (see Table I) the equality of 
Eq. holds within a few percent accuracy. This agreement is based on what all 
the CQMs have in common: the effective degrees of freedom of the three con- 
stituent quarks and the underlying SU(6) symmetry. 

On the other hand, the experimental value of the ratio is not reproduced by 
CQMs, at variance with the calculations based on phenomenological parton dis- 
tributions reported in Fig. 1. This means that the SU(6)-breaking mechanism 
contained in the phenomenological partonic distributions does not correspond to 
the SU(6) breaking mechanism implemented in the CQMs we have analyzed. 
The quark densities as given in Eqs. (|13],|14D are evaluated in the rest frame, as 
we are using non-relativistic wavefunctions in this paper. It is clear that in this 
way, relativistic boost effects are not included. Further work to quantify these 
relativistic boost effects is underway, even if we do not expect them to change in 
any important way our conclusion for the ratio of Eq. (|lD. 

To conclude, it seems that all CQMs are too strongly constrained by the pres- 
ence of the standard degrees of freedom corresponding to three constituent quarks. 
Therefore additional degrees of freedom should be introduced, in particular quark 
antiquark pairs and/or gluons and the discussed equation of Ref. [||, being sen- 
sitive to the SU(6)-breaking mechanism, will provide a useful tool for testing the 
new models. 
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Figure 1: Scale dependence of the rhs of Eq. (||) for various phenomenological forward 
parton distributions as indicated on the curves. Dotted curves : MRST parton distribu- 
tions (MRST98 NLO ||], MRSTOl NLO [||], MRSTOl NNLO |7|). Dashed curves : 



CTEQ parton distributions (CTEQ5M NLO |0|, CTEQ6M NLO Dashed-dotted 
curve : GRV98 NLO (MS) Also shown is the Ihs of Eq. (|), i.e. the experimental 
value for k^'/k" (constant solid curve). 
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Figure 2: Scale dependence of the proton momentum fraction Mf""' calculated with 
the Isgur-Karl model. 
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-1.0 


-1.0 


-1.0 


-0.9372 


R-ratio at = 0.5 GeV^ 


-1.0098 


-1.0030 


-0.9983 


-0.9881 


R-ratio at = 5.0 GcV^ 


-1.0098 


-1.0030 


-0.9983 


-0.9881 


R-ratio at = 10.0 GeV^ 


-1.0098 


-1.0030 


-0.9983 


-0.9881 



Table 1: Different CQM predictions for the R-ratio and for the K-ratio k^/k' 
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Figure 3: The R-ratio for the HCQM with isospin dependence compared with 
the K-ratio calculated with the same model and with the experimental value of 
the K-ratio . 
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